Singularity-Free Cylindrical Cosmological Model 
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Abstract 

A cylindrically symmetric perfect fluid spacetime with no curvature 
singularity is shown. The equation of state for the perfect fluid is that 
of a stiff fluid. The metric is diagonal and non-separable in comoving 
coordinates for the fluid. It is proven that the spacetime is geodesically 

■ complete and globally hyperbolic. 
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CT 1 - 1 Introduction 

u 

Due to the powerful singularity theorems (cfr. for instance pQ) it was widely 
believed that cosmological models were to have an initial singularity. However 
in 1990 Senovilla showed the first cosmological perfect fluid solution of the 
r> | Einstein equations with regular scalar curvature invariants. It corresponds to a 

■ cylindrically symmetric spacetime filled with an isotropic radiation perfect fluid. 
In [3j it was shown that this spacetime is indeed singularity-free. This is in full 
accordance with the singularity theorems, as it is to be expected, since although 
it fulfills the energy, generic and causality conditions, it does not have have any 
of the causally trapped sets required by the theorems (closed trapped surfaces, 
compact achronous sets without edge,...). This special solution was generalized 
in a subsequent paper where an Ansatz of separability of variables for diago- 
nal orthogonally transitive commuting G2 metrics was extensively explored. A 
thorough discussion of this Ansatz can be found in 0. This family is shown 
to be included in a wider class of separable cosmological models which com- 
prises FLRW universes [B]. Other properties of these solutions, such as their 
inflationary behaviour, their generalized Hubble law or the feasability of con- 
structing a realistic non-singular cosmological model, are studied therein. Other 
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non-singular cylindrically symmetric perfect fluid spacctimes have been found 
0, [El- Together with some solutions in these are the only non-diagonal 
ones that are known. Furthermore, in the family of solutions included in |1(J|. 
there are non-separable non-singular perfect fluid spacetimes, with and without 
symmetry axis 

In this paper we present a cylindrically symmetric stiff perfect fluid solution 
which is non-separable in comoving coordinates and has non-singular scalar 
curvature invariants. Due to these facts it is not included in the families quoted 
before. The physical properties of these metrics and their relation to other 
families of solutions are investigated. 



2 The metric 

The line element for a spacetime that admits an Abelian two-dimensional or- 
thogonally transitive group of isometries acting on spacelike surfaces can be cast 
in the form 

ds 2 = e K{t ^ (-dt 2 + dr 2 ) + e~ u ^ dz 2 + e u ^ r 2 d^ 2 , (1) 

where <f> and z are coordinates adapted to the commuting generators (in this 
chart, £ = £ = d z ) of the group. The remaining coordinates, {t,r}, have 
been chosen so that the line element induced in the subspaces orthogonal to the 
Killing orbits is isotropic in these coordinates. 

Following 11 , if we are to have a regular symmetry axis on the locus where 
A = <?(£, £) = 0, then we have to impose that 

ff(gradA,gradA) 

4A * L ' [2) 

on approaching the axis. 

The metric functions for the spacetime that we want to show are given in 
this set of coordinates by 

K(t,r) = r 4 + (a + (3)r 2 + 2t 2 (3 + At 2 1 r 2 (3) 

U(t,r) =/3(r 2 + 2t 2 ) (4) 

- oo < i, 2 < oo, < r < oo, O<0<27r (5) 

and are easily checked to satisfy the regularity condition on the set of events 
defined by r = 0. There are no further isometries than the ones that have 
been implemented from the beginning and therefore this spacetime is properly 
said to admit cylindrical symmetry. The only restriction that we impose on the 
parameters a and is that both of them are positive. 

The matter content in this spacetime is a stiff perfect fluid, whose density 
and pressure, 

f i = p = ae- K ^ r \ (6) 
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are regular everywhere in this chart. When a equals zero we have a vacuum 
solution, from which the stiff fluid metric can be recovered by means of the 
Wainwright, Ince, Marshman algorithm However this solution has not 

been obtained by this method, although it could be used as it is done in |5] to 
generate further models with the same vacuum metric. 
The velocity of the fluid takes the form, 

u = e-i K ^d t , (7) 

and therefore the coordinates are comoving. From the expression for the line 
element Q it is clear that it cannot be rendered separable in comoving coordi- 
nates. 

The acceleration of the fluid has projection only on the radial direction, 

a = r {[3 2 r 2 + a + (3 + 4 f3 2 1 2 ) d r , (8) 

due to the orthogonal transitivity requirement and that the fact that the velocity 
is orthogonal to the orbits of the group of isometries. For the same reason the 
vorticity is zero. 

In an orthonormal coframe {8°, 9 1 , 8 2 , 8 3 }, where the four independent dif- 
ferential one-forms take the form, 



go = e * d^e^^dr, 8 2 = <T* U ^ dz, 9 3 = e^'^r cfy, 

(9) 

the shear tensor constructed with the derivatives of the velocity u can be written 

as, 

a = -f3tc-^ K ^ {( 1 + 2 f3r 2 ) 8 1 g> 8 1 - ( 2 + f3r 2 ) 8 2 ® 8 2 + (1 - f3r 2 ) 8 3 ® 8 3 } 

(10) 

The expression for the expansion, Q, of the cosmological fluid, 

e = 2/3t(l + 2/3r 2 )e-* if( *' r) , (11) 
allows us to calculate the deceleration parameter, q, for this universe, 

Q = 2 ~ 2 /3i 2 (l + 2/3r2)' (12) 
from the action of the vector field u on the inverse of <d, ^Sj 



it 



S)- 5 (! + «)• (13) 



Since q is positive in the time span (— ti„/, U n f), 



tmf ~ V4/3 (1 + 2/3 r 2 )' (W) 
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we are led to conclude that this spacetime has an inflationary epoch, which is 
longer the closer the observer is to the symmetry axis. Furthermore the scale 
factor of the universe can be defined ^3] as a solution of the differential equation, 

u(R)=®R, (15) 

which in our case can be solved, introducing an arbitrary function C of the 
spacelike coordinates, 

R( t, r, z,<t> ) = C(r, z, <£) e^ * <■ 1+2 r " > . (16) 

3 Regularity of the metric 

Since the scalar invariants that can be formed with the metric and the Riemann 
curvature are polynomials of the density and the pressure of the fluid and the 
components of the Weyl tensor, we shall calculate the latter ones in order to 
show that the curvature scalars are not singular. In the null tetrad that can be 
naturally constructed with the one- forms of ©, the components of the Weyl 



tensor can be shown to be, 

^ = \{fx{t,r) + h{t,r))e- K ^ (17) 

*i=0 (18) 

* 2 = \ (3 (3 2 r 2 + 3 - a - 12 2 t 2 ) e~ K ^ (19) 

* 3 = (20) 

*4 = l(fi(t,r)-f 2 (t,r))e- K «> r \ (21) 

/i^.r) = -30 + 20 3 r 4 + 30 2 r 2 + a(l + 20r 2 ) + 2A0 3 r 2 t 2 + 120 2 t 2 (22) 

/ 2 (t,r) = 12 P 3 r 3 t + 12 (3 2 rt + 16 (3 3 1 3 r + 4a (3tr. (23) 



From the expressions for the components of the Weyl tensor and the density 
and the pressure of the fluid it follows that they vanish when either t or r tend to 
infinity. Hence the spacetime is flat for large values of the time and radial coor- 
dinate. It is also easy to check that the Weyl components are regular everywhere 
and so are the pressure and the density. Therefore all the curvature invariants 
are regular. The spacetime has a low matter content for large negative values of 
the time coordinate and undergoes a contracting epoch until t reaches the zero 
value. For positive values of t this universe is expanding. The change from con- 
traction to expansion without developing a curvature singularity appears also 
in other non-singular cosmological models 0, [H], [I]. 
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The gradient of the comoving time coordinate t is always negative. It is 
therefore a cosmic time £Q and the spacetime is causally stable. In particular 
this implies weaker causality conditions suchs as the chronology condition. 

The strong and dominant energy conditions are satisfied since the density 
of the fluid is positive everywhere and the equation of state corresponds to a 
stiff fluid. Moreover, the energy-momentum tensor does not vanish anywhere. 
This means that for every non-spacelike vector W the contraction with the Ricci 
tensor R(W, W) is greater than zero and it implies that the generic condition is 
satisfied as well [Tl]. 

There is a static limit for the metric which amounts to take the parameter 
[3 equal to zero. The resulting metric, 

ds 2 = c ar2 (-dt 2 + dr 2 ) + dz 2 + r 2 dcj> 2 , (24) 
can be seen to be also the static limit of 0. 

4 Geodesic completeness 

In order to determine whether the spacetime is geodesically complete, we have 
to study the equations for the causal geodesies, 

W v v = 0, v — id t + r d r + z d z + (pd,p, (25) 

where v is the tangent vector along the geodesic and the dot stands for the 
derivative with respect to the affine parameter r. We shall introduce a constant 
of motion 8, which takes the zero value for null geodesies and one for timelike 
geodesies, that is, 

g(v,v) = -S. (26) 

The existence of isometries simplifies the problem, since two new indepen- 
dent constants of motion arise, 

L = g(v,d 4> )=r 2 e u ^4>, Z = g(v, d x ) = e~ u ^ z, (27) 

as a consequence of the geodesic and the Killing equation. 
Since 121)1 can be cast in the form, 

i 2 - r 2 = e -K(*,r) {5 + Z 2 e u ^ + L 2 e~ u ^ r~ 2 }, (28) 

it is convenient to parametrize t and r with hyperbolic functions of a new vari- 
able, £, as it is done in 0, in order to lower the order of the geodesic equations. 
The final system of equations for future-directed geodesies (past-directed ones 
can be handled in a similar way) is first order, 



-e 



ijC(t,r) J r cosH (Z 2 A(t, r) + L 2 B(t, r) + 6 C(t, rj) 

{ 7W7) + (29) 
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t sinhg (Z 2 D(t,r)+L 2 E(t,r) + S F(t,r)) 1 

+ VW^) J 

i = e K & r y 2 G(t,r) cosh£ (30) 
f = e K ^ r ^ 2 G(t,r) sinh£ (31) 
i = c 6 ^^ Z (32) 

^ = e- a ^r- 2 i, (33) 

and the last two equations are just quadratures that can be integrated after 
solving the first three equations. We have introduced seven functions of t and 

r, 

A(t, r) = (/3 2 r 2 +A/3 2 t 2 +a + 2 0) e u ^ (34) 

S(*,r) = (^-l + ^±^)e-^) (35) 

C(t,r) = (3 2 r 2 + a + P + Ap 2 t 2 (36) 
D(t,r) =4/3(1 +(3r 2 )e u ^ r) (37) 
E(t,r) =A/3 2 e- u ^ (38) 

F(t,r) = 2/3 (1 + 2/3 r 2 ) (39) 

f2.-[/(t,r) 

G(i,r) =5 + zV / ^ + ----^ . (40) 

The only geodesies that fall out of this scheme are radial null geodesies, 
S = L = Z = 0. For them the analysis is fairly simple, since t = \r\ and the 
only equation that is left for integration is, 

r + 2r(/3 2 r 2 + a + /3 + 4/3 2 i 2 ) f 2 + 4 /3< (1 + 2 /3r 2 ) if = 0, (41) 
which can be seen to have a first integral, h, 

h = f c 2 fit 2 (1+2 Pr 2 ) + (a+fl)r 2 + ±f3 2 r\ ^ 

The completeness of these geodesies is obvious, since r is confined finite in 
the interval [— \h\, \h\]. 

Going back to the generic equations [2^1 we shall explore the appearance of 
divergences that may occur when the radius, r is either too small or too large 
or when the time coordinate, t, is too large. 

If the radius, r, grows too large, so that r 4 > /3~ 2 , r > (sinh^ > 0), then £ 
is negative when the time coordinate, t, is positive. Therefore £ decreases and r 
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does not diverge for finite affine parameter since r cannot grow arbitrarily large. 
If t is negative, since the time coordinate grows at least as fast as r, t becomes 
positive before r diverges. 

The terms with negative powers of r do not yield singularities when the 
geodesies approach the axis (f < 0, sinh£ < 0). The reason is that £ becomes 
positive when the geodesic is close to r = and prevents the radial coordinate 
from decreasing too quickly. This happens because the only negative terms in 
the expression of £ are either bounded or overcome by the term in r~ 4 in B(t, r) 
for decreasing radial coordinate. 

Also if t grew too fast, this would mean that r would grow or decrease quickly 
and the previous reasonings would prevent divergences from appearing. 

Hence we are led to conclude that every causal geodesic is complete and 
therefore the spacetime is geodesically complete. 

From this analysis it follows that every null geodesic intersects once and only 
once every timelike hypersurface t — const. Therefore these hypersurfaces are 
Cauchy surfaces and the spacetime is globally hyperbolic |15) . 

5 Discussion 

It has been shown that the spacetime is singularity-free and globally hyperbolic. 
Since energy, generic and causal conditions are fulfilled, it is clear that this 
spacetime does not have causally trapped sets £Q, The gradient of the 

transitivity surface element is spacelike as in all known non-singular models. 

The question of whether there is an open set of non-singular cosmological 
models in a reasonable topology remains open. However the fact that both 
this spacetime and the one shown in [7j have the same static limit suggest that 
more solutions could be found by suitably adding new parameters to the static 
solution. 
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